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ABSTRACT 

In order to establish a model which would more accurately predict 
discharge characteristics of a cell at any arbitrary current, and which 
would better fit the family of discharge curves, testing was done based 
on the Shepherd equation using a transparent lead-acid battery. 

A set of discharge curves at .6, 1.5, 3.6, and 5-^A, preceded by a 
3.6 a preconditioning discharge to remove the ’’memory effect," was obtained. 
The battery was cycled continuously through deep discharges in an attempt 
to determine secular changes during the life of the battery. 

Several modifications to the Shepherd equation were proposed in 
order to obtain a better fit to discharge data. These included modeling 
the capacity as a function of current as Q=Cl(l“ n ), Peukert's equation; 
modeling the ohmic polarization as a linear function of accumulated charge; 
and modeling the diffusion polarization as a function of accumulated charge 
alone and not current. 

Conclusions were drawn on the suitability of Shepherd’s equation. 
The inclusion of the changed ohmic polarization term gave a good fit for 
individual curves which also suggested the change in the diffusion polari- 
zation term. Since a small current density approximation was used in 
Shepherd’s equation, developing a model based on the original terms was 
recommended. Also recommended were the modification of a model to predict 
characteristics during a cycle with various discharge currents and further 
life cycle testing to determine the secular changes in the battery. 
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CHAPTER 1 



INTRODUCTION 

1.1 BACKGROUND 

Predicting the transient and long-term characteristics 
of a cell under arbitrary charging and discharging histories 
is a problem confronting a system designer which could be 
overcome if there were an electrochemical model of the 
lead-acid cell. With this model, it would theoretically 
become possible to determine the optimum scheme for getting 
energy into and out of a battery. 

To what degree is the efficiency of a battery affected 
by different types of discharge? Should the designer select 
a pulse-type discharge over a steady discharge? What effect 
does this have on the system as a whole? To what degree does 
the age of the battery affect the efficiency? Low current or 
high? What exactly is the effect of temperature and current 
on the charging and discharging history of the arbitrary 
battery? Will current models handle the problem of the age 
of the battery? Does the model provide a fit merely to a 
single discharge curve or for the family of curves? Does it 
account for the condition of the battery? What effect has 
gassing on the capacity for full charge or discharge? 
Questions such as these pose a great problem for the designer 
of a model. 

Naturally, there are numerous other criteria which might 
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affect the optimum charge/discharge history of a battery, so 
it becomes very important to select a battery type to test 
which possesses characteristics which are relatively similar 
to other cells. 

The lead-acid cell was chosen for this investigation for 
several reasons. Lead-acid batteries are the most widely 
used batteries in the world, and for that reason, they 
represent the greatest commercial interest. Portable lead- 
acid batteries have automotive use, use in motorcycles, 
railroads, all types of military vehicles on the land, sea, 
and air, and additional uses not related to vehicles. 

Additionally, there is a wealth of specific information 
concerning lead-acid batteries in literature. Volume after 
volume has been written exploring the numerous characteristics 
of this type of battery in great detail. Chapters upon 
chapters have been written about the history, electrochemistry, 
applications and performance characteristics of lead-acid 
batteries. Detailed studies have been done on performance and 
structural aspects, as well as on self-discharge effects. 

Because the author chose as the aim of his thesis, the 
establishment of an electrochemical model of the cell, based 
on Shepherd's equation (1965), it was necessary to chose a 
battery type for the experimentation which was rechargable, 
lightweight, and portable, in addition to having similar 
characteristics of other types of cells. The lead-acid 
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battery, in addition to being rechargable, is available in a 
clear plastic case, thus affording the investigator a visual 
means to determine some battery conditions which are not 
readily available empirically in a battery with an opaque case 
Also, it was felt that the characteristics of some of 
the other types of cells are similar to those of the lead-acid 
cell and that a model for the one may be applicable to the 
others. 

1.2 HISTORY OF LEAD-ACID BATTERIES 

The original incentive toward the development of an 
electrochemical source of current was the "galvanic cell" 
which was based on Galvani's "frog leg" experiment, circa 
1789. 

Volta, after 1792, developed several useful current 
sources, such as the Volta pile. He also formulated the 
first electrochemical theory on the subject. 

Ritter developed a cell in about 1803, Ritter's pile, 
which, when disconnected from an external source of current, 
would yield a current in the opposite direction. This was the 
basis for the secondary cell or storage battery. 

In 1860, Plante demonstrated his cell in which two lead 
sheets separated by rubber were rolled into a spiral and 
placed in a 10% sulfuric acid solution. When charged, the 
amount of energy depended upon the amount of lead oxide formed 
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By 1880, W.V. Siemens had invented the dynamo, which 
was then installed in central power stations. Batteries were 
used for load leveling purposes, being charged with energy 
produced at night. 

Around 1881, Faure coated the lead foil with lead oxide 
for the purpose of accelerating Plante's forming process. 

He used various lead current collector grids. 

1883 brought us the Plante plate which was actually 
designed by Tudor. In this design, he pasted lead oxide onto 
cast ribbed pure lead plates. Grid plates, plates cast from 
antimonial alloys were also developed at about the same time. 

From 1890 to 1920, wood was used to separate the plates 
and then rubber. Around 1927, porous ebonite and microporous 
rubber was used as separators. 

Around 1965, batteries were designed to be relatively 
maintenance-free for use in portable items, and, with the 
advent of plastic battery cases and lightweight construction, 
automation and mass-production became feasible. 

The 1970s gave us new applications of the lead-acid 
battery, with increasing interest in Traction batteries coming 
to the forefront. Water-activated and maintenance-free 
automobile batteries were designed. Kordesch, (JL977) 

These historical developments, some important, some 
merely historically interesting, represent the basis for the 
lead-acid battery as it is known today. 



1.3 THE NEED FOR A MODEL 



A system designer needs an accurate electrochemical 
battery model in order to produce a design which makes optimum 
use of these batteries. With the use of computer-aided 
design, a functional description of various battery character- 
istics would greatly simplify the present cumbersome task of 
establishing an algorithm for interpolation of battery 
characteristics from its family of discharge curves. An 
algorithm for predicting the response to an arbitrary (non- 
constant current) discharge or charge would be, at the very 
least, extremely difficult, if not prone to significant 
inaccuracies . 

Several empirical models have been developed in the past 
to show the dependence of capacity on discharge current I 
and discharge time t {Schroder (1894), Liebenow (1897), 

Peukert (1897), Rabl (1936), Davytan (1946), and Shepherd 
(1965)}. Shepherd's (1965) model (discussed in detail in 
Chapter II) gives a relationship between cell voltage on dis- 
charge as a function of current density and time based on the 
theoretical laws of electrochemistry. Shepherd's equation 
(1965) serves as the basis for this investigation into estab- 
lishing a State Model. 

Very early in this investigation, it was found that the 
equation provided a satisfactory fit to a single discharge 
curve (complete discharge at constant current) but a poor fit 



for the family of discharge curves. A typical family of 
discharge curves is shown in figure 1-1. Also, Shepherd 
did not include any term which would account for secular life 
changes in the cell. 




Fig. 1-1. Discharge Curves for a Typical 100-Ah (20-h rate) 
SLI Lead-Acid Battery. 

The emphasis of this thesis was to establish a model 
which would more accurately predict the discharge characteris 
tics of a cell at any arbitrary current, thereby providing a 
better fit for the family of discharge curves. Basing his 
investigation into a State Model on the Shepherd equation 
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(1965) , the author has chosen to call his model the "Modified 
Shepherd Equation." It was also the desire of the author 
to experimentally determine the secular life changes of the 
lead-acid cell and to include a model of these changes in the 
Modified Shepherd Equation. 

The problem of needing an adequate model was made clear 
to the author while working on a Diesel submarine computer 
synthesis model at the Naval Ships Engineering Center in 
Washington, D.C. A great deal of time was invested in estab- 
lishing a model of the lead-acid cell which would adequately 
predict the capacity of the cell at arbitrary discharge cur- 
rents. This was accomplished simply as a curve fit to the 
capacities determined from the family of discharge curves. 

Due to the complexity of the problem, no attempt was made 
to use the discharge characteristics in that model. If a model 
such as the one proposed in this thesis had been available, not 
only would a significant number of man-hours have been saved, 
but also another degree of flexibility would have been allowed 
in the synthesis model. 

Also, the problem of the design of an electric vehicle is 
essentially identical in this respect. Moreover, the battery 
manufacturer, by supplying only a few coefficients for a 
battery model would be providing the characteristics of his 
product in a very concise form which would require no inter- 
polation for the user's specific application. 
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1 . 4 BATTERY USAGE 

Since the earliest "galvanic cell", the lead-acid storage 
battery has been in demand for numerous purposes. Marine 
uses alone are diverse. 

Marine batteries, which are rugged and durable, have 
given long service to ships of the Navy and Merchant Marine. 
These batteries function under extremely adverse conditions 
and are subject to continual vibration. 

These batteries have been used in all types of ships from 
the largest aircraft carrier to the smallest pleasure craft to 
power radios, lights, pumps, hoists, and compressors. They 
are used in torpedoes and to fire guns. They are used in the 
kitchens of ships to provide power to stoves and refrigerators. 

Batteries power alarms and emergency lighting and even 
pump bilges and hoist anchors. Many are used in ignition 
circuits of internal-combustion engines of small craft. They 
provide emergency power for submerged Nuclear submarines and 
the propulsion for submerged Diesel submarines. 

Marine batteries are used for gyrocompass emergency power 
and for telephones, and emergency generator start-up/ or 
field flashing. 

For many years , scientists and engineers have studied the 
possibility of converting from gasoline powered motor vehicles 
to electrically powered motor vehicles, from hybrids to all- 
electric delivery vans. Much interest was in all-electric 



automobiles. It has been found that there are no savings in 
energy if batteries simply replace the gasoline engine. 
Electrically powered vehicles would offer a significant 
pollution advantage over gasoline and Diesel engines. 
Batteries could be used to provide additional energy require- 
ments in hybrid gasoline/electric vehicles for load leveling 
purposes or for applications where the possible recovery of 
energy through regenerative breaking is significant. 

(Kordesch 1977) 

Batteries have been classified according to usage, 
according to Bode (1977) : 

1. SLI Batteries (Starting, Lighting and Ignition) 
These batteries are used to start an internal- 

combustion engine and to provide power to the 
electrical system when the engine is not running. 

The largest number of storage batteries are for 
use in the automobile and are used for this purpose 
more than for any other . 

2. Stationary Batteries 

These generally provide emergency power when 
the primary source of power is not available. 

These include such applications as communications 
systems, electric utilities, computer systems, 
emergency lighting, and railways to provide peak 
loads, emergency power or smoothing of ac power. 
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3. Motive Power (Traction) Batteries 

Lead storage batteries provide power for the 
propulsion of devices such as electric lift trucks, 
mining equipment, street delivery vehicles, and 
other types of material handling. 

4. Special Purpose Batteries 

Some examples of batteries in this category 
are aircraft, submarine, special military and small 
sealed batteries for consumer applications. Some of 
these uses have already been discussed. 
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1.5 BATTERY LIFE 

The life of a battery can be defined as the period of 
time during which it has the ability to perform within certain 
specified limits. In other words, the end of life occurs 
when the battery can no longer perform within these limits. 
These limits actually depend upon the application and the re- 
quirements of the system they serve. The end of life of an 
automobile battery occurs when it no longer is able to provide 
sufficient cranking power to start the engine. For a battery 
powered vehicle, it may be when the vehicle will no longer 
travel a certain number of miles at a specified speed; in 
traction batteries, the decrease of capacity below a certain 
value may be the criterion on which the end of life decision is 
based. Generally, however, systems are designed to meet all 
specified requirements as long as the battery is able to 
deliver 80% or more of its new capacity. 

Bode (1977) states that the end of battery life is best 
evaluated in the laboratory using fixed standards and speci- 
fied test conditions. The ideal requirement is that the life 
test correspond to the load conditions of the battery in the 
field, however, as in the case of automobile batteries, ideal 
testing in the laboratory is impossible because it is im- 
possible to duplicate such diverse conditions as automobile 
type, driving methods, maintenance and environment. 

Figure 1-2 is a plot of the change in capacity during the 



14 



lifetime of a lead-acid battery. Very little information is 
available which provides details on the change of characteris- 
tics of a cell during its lifetime. It is therefore proposed 
that this data be obtained by experimental means. 




Figure 1-2. Typical Cell Capacity During the Lifetime of 
a Lead-Acid Battery. KORDESCH (1977) 

The discharge capacity of a battery has been found to be 
somewhat dependent upon the rate of previous discharges. Bode 
(1977) calls this the "memory effect." Vinal (1955) states 
that "the capacity is lower if the discharge has been preceded 
by a discharge at a higher rate, and it is higher if preceded 
by a discharge at a lower rate." Therefore, any experimental 
procedures should take the "memory effect" into consideration. 
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1.6 REASONS FOR DECREASE IN LIFE EXPECTANCY OR FAILURE 

According to Vinal (1955) , storage batteries large and 
small and in widely diverse services are subject to conditions 
that may lead to excellent service or through accident or 
abuse to something less than satisfactory. Under his heading 
"Sources of Trouble," he writes about faulty conditions 
which if detected could be corrected. Detection and correction 
of faulty conditions might prolong the satisfactory battery 
life into years. 

He also states that as a result of studies done in 1947 
on thousands of worn out batteries which had been used in 
normal service, the reasons they wear out became well known . 

In comprehensive data relating primarily to automobile 
batteries, grid corrosion of the positive plates accounted for 
42% of the total failures. Naturally, with the use of new 
alloys, batteries developed since that time should have less 
positive plate corrosion. 

Bode (1977) indicates that destruction of the positive 
grids occur when lead sulfate forms a heavy layer on the grid 
to protect it from being dissolved anodically. One rarely 
encounters corrosion of the negative grid. Hydrogen corrosion 
occurs in the presence of organic substances and at higher 
operating temperatures. The result of corrosion is loss of 
conductivity and grid breakage. 
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Vinal (1955) indicated that 7 to 10% had less than total 
life span due to shedding of positive active material, 
cracked partitions and leaking cases accounted for nearly one- 
third, and buckled plates and short-circuited separators pro- 
duced failure in 10 to 12% of the batteries. About 10% of the 
failures were caused by negative plate failure. Sulfation 
which causes buckled plates and negative failures, was a cause 
of battery failure in less than 20% of the batteries. 

Overcharging loosens active material, corrodes positive 
grids, causes excessive gassing and produces sludging of a 
fine brown sediment. It may cause raised temperatures contri- 
buting to battery destruction with damage to plates and separ- 
ators, buckling of plates and needless loss of water. 

Bode (1977) states that in sulfation, . "hard" lead 
sulfate forms if chargeable lead sulfate undergoes long or high 
temperature aging process, (recrystallization). Sulfation 
produces loss of capacity and is caused by long standing in a 
discharge condition, too high acid concentration, increased 
self-discharge , continuous operation between 45° and 50°C, 
and prolonged undercharging. 

He also indicates that sludging of the positive mass 
occurs most at the end of a charge and at the beginning of a 
discharge. Heavy gassing facilitates sludging. There is a 
greater amount of sludging in discharges with large current 
density. The type of discharge influences sludging. If acid 
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concentration is drastically reduced, temperature increases or 
current density is diminished, there is less sludging and 
the battery life is increased. 

Vinal (1955) concludes that occasional overcharging is 
beneficial, but habitual overcharging decreases useful battery 
life. 

He indicates that undercharging gradually runs down cells 
producing diminishing values of specific gravity. Plates be- 
come lighter in color and a sludge, a fine white powder which 
is predominantly lead sulfate, is produced. One or more cells 
could become exhausted. He concludes that the remedy is to 
charge until all cells are in normal condition. This is refer- 
red to as an equalizing charge. Undercharging is the most 
common cause of buckling of plates because lead sulfate occu- 
pies more space than the original material producing excessive 
strain on the plates. 

Failure of the charging system and of the battery ulti- 
mately is blamed on the battery. When an automobile reaches 
speeds greater than or equal to 12 to 15 mph the battery begins 
to charge. If the generator doesn't come on at the proper 
speed the battery won't receive its required amount of charge. 
Also, the system must regulate itself properly so that the 
battery will not receive massive amounts of current when the 
automobile reaches highway speeds. Conditions in the charging 



28 



system which drain the battery and cause it to remain dis- 
charged are grounded circuitry and partial grounds. 

(Vinal , 1955) 

Other possible causes of shortened battery life include 
corroded terminals, short circuits, separator failures, 
worn-out plates, electrolyte below the tops of the plates, 
freezing, when the battery is in a discharge condition, 
impurities in the solution, spalling of the active material, 
growth of positive plates, shrinkage of negative plates, 
leading, and battery explosions. The interested reader is 
directed to the classic text of Vinal (1955) and the excellent 
texts of Bode (1977) and Kordesch (1977) . 
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CHAPTER 2 

THEORETICAL ANALYSIS 

2.1 ELECTROCHEMISTRY OF THE LEAD-ACID BATTERY 

The theory of operation of the lead-acid battery is 
based on the "double sulfate theory" proposed by Gladstone 
and Tribe, (1881-1883) . 

Pb0 2 + Pb + 2H 2 S0 4 dlschar 9 a 2PbS0 4 + 2H 2 0 (2.1) 

charge 

The principle was so named because the discharge product 
of both the lead anode electrochemical oxidation reaction and 
the lead dioxide cathode electrochemical reduction reaction is 
lead sulfate. 

This early equation reveals the stoichiometric relation- 
ship of reactants and products, as well as the role played by 
the sulfuric acid. 

A more complete representation of the actual electro- 
chemical reactions is: 

ANODE: Pb + HSO^ dlSCharg 1 PbS0 4 + H + + e“ (2.2) 

charge 

CATHODE: Pb0 2 + 3H + + HSO^ + 2e~ - » PbS0 4 + 2H 2 0 (2.3) 

OVERALL: Pb + Pb0 2 + 2H + + 2HSoZ - » 2PbS0 4 + 2H 2 0 (2.4) 









' 
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The electrolytic dissociation of sulfuric acid into 

H + and HSO 4 is closer to the true picture in the cell. The 

+ — 2 

dissociation into 2H and SO4 occurs in solutions much more 
dilute than those found in a lead-acid cell. 

During overcharge, the electrolysis of water occurs, 
yielding hydrogen and exygen. 

Based on the free-energy change of 

AG = -13.95 kcal/mole 

for the anodic reaction 

(E q = 0.303 V versus H 2 /H + ) 



and 



AG = -74.99 kcal/mole 
for the cathodic reaction 

CE q = -1.627 V versus H 2 /H + ), 
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the overall reaction values are: 

AG = -88.94 kcal/mole 
CE fl = 1.93 V) 

specific energy = 161 Wh/kg 
specific capacity = 83.5 Ah/kg 

Bode (1977) . 

In practice, significantly lower specific energies are 
exhibited, but the theoretical values provide an upper limit to 
possible future developments. For a complete discussion of 
electro-chemical theory, the reader is directed to the 
excellent text by Bode (1977) , and the classic text in the area 
by Vinal (1955) . 



32 



2.2 THE SHEPHERD EQUATION 

Shepherd (1965) has proposed a static model for cell 
discharge characteristics which uses simple analytical 
expressions of chemicophysical origin and describes the cell 
potential E as follows : 



E(q,i) = E 0 - V a (q) - V 0 (i) - V d (i,q) 

where E 0 = open circuit cell potential 

V (.q) = A exp (- - ) = activation polarization 

Vo (.i) = Roi = ohmic polarization 



V d Ci,q) = K ( 



Q 

Q-q 



i = diffusion (concentration) 
polarization 



i = the discharge current 
q = it = the accumulated charge 
Q = the cell capacity 
K = the mean polarization resistance 
Ro = the internal cell resistance 



Shepherd used "it" in place of "q" in his equation 
and R 0 corresponds to his N. The activation polarization 
term is the single term which does not have any particular 
theoretical basis. Its inclusion accounts for the rapid 
drop in potential at the beginning of discharge. Figure 2-1 
shows an exagerated view of this phenomenon. The time 
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duration of this initial drop was found to be in order of se- 
conds. Consequently, the description of cell discharge should 
not be significantly affected by simply ignoring the term; this 
was the approach used by Shepherd, also. 

Another phenomenon shown in fig. 2-1, not treated by Shep- 
herd, is known as the "voltage dip" or "coup de fouet" where 
initially the voltage drops below the "steady state" value and 
then rises to this level within a few minutes. 




Figure 2-1. Cell Voltage versus Time for a Constant 
Current Discharge 

The ohmic polarization term accounts for the ohmic resis- 
tance through both the metal and the solution. This has been 
modeled as a constant, however, it is felt that this assump- 
tion may be erroneous. 
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The specific electrical resistance of lead oxide is on 
the order of 10“ 6 £2/m and metallic lead, on the order of 
10“ 7 £2/ m. The electrolyte has a range of specific resistance 
of about 125 £2/m when nearly fully charged to about 17 0 £2/m 
when discharged. The product of the two reactions, lead 
sulfate, is highly resistive, specific resistance on the 
order of 10 10 £2/ m. Based on these values, it would be expect- 
ed that the ohmic resistance increase as the reaction proceeds 
and becomes fairly significant as it approaches completion. 

The diffusion (concentration) polarization term appears 
to have been made up from a combination of what Bode (1977) 
calls the charge transfer overvoltage, n. , and the concen- 
tration overvoltage, n . Following Bode's development, 
the terminal voltage of the cell is 



E 



E 

rev 



+ 




(n t + Vc ” IR 



( 2 . 6 ) 



R rev ; '' s °P en c i rcu it voltage and IR includes the 

resistances of the electrolyte and separator. The subscripts 
A and C refer to the anode and cathode overvoltages, respec- 
tively. All of the overvoltages decrease the cell voltage on 
discharge, and overvoltages of the electrodes are related to 
the current densities, i. Assuming that the electrodes have 
equal surface area A a = A p = A, then for both electrodes the 



35 



following expressions apply: 

(2.7) 

( 2 . 8 ) 

For small current densities, these become 



n t = 2 ‘ 3 RT log i 0 (1 -0 ) - b log i 



and 



zF 

2.3 RT 



n c zF 



log (1 - 



) . 



1 *. = 



RT 



zFi o (1 -0) 



and 



ri = 



RT 



zF 



(2.9) 



( 2 . 10 ) 



For this case, one of the terms in (2.6) may be neglect- 
ed if they differ by an order of magnitude. The terms in the 
equations above are: 

i 0 = the apparent exchange density 
z = the number of electrons transferred in the 
determining step 

F = the Faraday (measure of capacity) 

R = the gas constant 
T = the temperature 

i^ = the diffusion limiting current 

and 

0 = the degree of coverage. (only the uncovered portion 
(1-0) of the electrode can react) 
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Replacing the constant terms in n t by a single constant k 
yields 

n. = (2.11) 

( 1 - 0 ) 

If the uncovered portion of the electrode is related to 
the cell capacity Q and accumulated charge q where q < Q 
with q + Q at complete discharge, equation (2.11) becomes 

_ _ K _ v Q . (2.12) 

r|, - - a 

(1-0) (Q-q) 

This term is similar to Shepherd's diffusion polariza- 
tion but with the striking difference that it is not a 
function of current. If one were to fit Shepherd's equation 
to a single discharge curve, this difference would not be 
apparent because the product of K and i in equation (2.5) 
would be equal to K in equation (2.12) . 

Prior to the discovery of the above relationship, the 
author, while attempting to make a fit for a family of dis- 
charge curves to Shepherd's equation, produced fits of each 
curve individually. In the process of analyzing these results 
a rather surprising discovery was made in that the product 
of K and i was nearly constant. 

Another item in Shepherd ' s equation which should be 
addressed, is the cell capacity, Q. The value for Q should be 
close to the capacity of the cell at the particular discharge 
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current. In fact, it needs to be slightly larger than the 
actual capacity to avoid predicting infinite overvoltages, 
i.e. Q > q for all values of q at the particular discharge 
current. 

The most widely used model for all capacity is Peukert's 
(1897) equation, 



In this equation n and C are constants which can be determined 
by using the discharge rate and time for complete discharge 
at two different currents. Since the cell capacity in terms 
of Q is the desired quantity, the substitution of t = Q/I gives 



(2.13) 




(2.14) 



The values of n and C can be calculated as follows: 



Q 



Qi 



2 




(n-1) log Ij = log C - log Q 



(n-l)log I 2 = log C - log Q 2 



log Q 2 - log Qj 



n 



+ 1 



log Ij - log I 2 



3S 



The value for n is in the range of from 1 to 2 . For 
diminishingly small current densities, n approaches 1 and 
for high current densities, n approaches 2. Typically, the 
value of n falls v/ithin about 1.3 to 1.4 This equation is 
not based on any particular physical laws, but serves as a 
good empirical fit to discharge capacities. 

Baikie, Gillibrand, and Peters (1972) expanded on this 
relationship to include temperature effects in the range of 
-18°C to +16°C; for current densities between 21 and 
165mA/cm 2 it has the form 

K 0 (1 + ct6) 

Q = (2.15) 

l ( n "0 

Where Q is the effective plate capacity, and 0 is the temp- 
erature in °C. The values they established were 

n = 1.4 for both plates 
K =0.32 and 0.24 min, 

and a = 0.021 and 0.015 min/°C 

for the negative and positive plates, respectively. 

Even if these equations provide a perfect correlation to 
the capacity of a cell, they do so only at one instant and 
do not predict changes which occur throughout the life of the 
battery. If a state model is to be established for the electro- 
chemical cell, then these effects must also be included. 
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2.3 METHOD OF DETERMINING COEFFICIENTS FOR THE MODEL 

An iterative approach using a digital computer was used to 
determine the coefficients for the model. This curve fitting 
method basically used the following steps: 

1. Choose initial values 

2. Vary each value independently 

n 

Trial Value = Value * (Factor) 
where n = -1,0,1 and 0 Factor 1 

3. Choose the set of trial values which gives the 
smallest difference between the data and model. 

Steps 2 and 3 were repeated until no further improvement 
was required by the user. When the "best" set of values was 
found, then these trial value were used in the next iteration 
as the values. When a particular factor failed to make any 
improvement, the difference between the factor and 1 was reduc- 
ed by 1/2 and a new iteration was performed. 

For runs of the Shepherd equation, initial values were 
chosen such that: 

E =. open circuit voltage 

Q = cell capacity plus a small amount to prevent 
division by zero 

Ki + Ri = E s - discharge voltage at t = 0. 
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This method finds a minimum value for the correlation, 
but the only assurance is that the set of values corresponds 
to a relative minimum. 
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CHAPTER 3 

EXPERIMENTAL PROCEDURE 
3.1 BATTERY CHARACTERISTICS 

In order to establish a model which would more accurately 
predict the discharge characteristics of a cell at any arbi- 
trary current and throughout its life, testing was performed 
using a lead-acid battery. 

The type of battery used was a low capacity motorcycle 
type battery in a clear case. The low capacity was used so 
that the equipment could be chosen requiring smaller current 
carrying capacities. The advantage of the clear case over an 
opaque case was that the author could more readily observe 
conditions of . the battery, including the amount of gassing 
during charging, the level of the electrolyte, the changes, if 
any, in the condition of the plates and the accumulation of 
any sludge, or fine-grained sediment, in the bottom of the bat- 
tery case without destroying the battery in order to make these 
observations . 

Table 3-1 lists the battery specifications and overall 
dimensions. Figure 3-1 shows the internal configuration of a 
cell. Each cell was configured with anode (-) plates on the 
ends; As can be seen, in the figure, there are four anode (-) 
plates and three cathode (+) plates with 6 separators in the 
spaces between the plates. 
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anode 

plate 



separator 



cathode 

plate 



28.2 mm 



- Anode Plates 




Figure 3-1. Internal Configuration of a Cell. 



Manufactured for: 



Battery Type: 

JIS Designation: 

Capacity: 

Case: 

External Case Dimensions: 



Wisco Division 
ESB Incorporated 
1222 18th’ Street 
Facine, VJI 53403 

JA-11 6V 3 Cell 

6N6-3B 

6Ah @ 10-h rate 

Clear Plastic Bottom 
with Black Plastic Top 

L = 10.7cm 
W = 5.5cm 
H = 11.0cm 



Table 3-1. Battery Specifications- and Dimensions 



4 3 




Figure 3-2. Side View of Anode Plate #1 as Seen from the 
End of the Battery. 



Table 3-2 lists the measurements on cell #1 to de- 
termine the thicknesses of the plates and separators. These 
measurements were made by placing a rule against the side of 
the case. Some degree of uncertainty in the measurement of the 
plate thicknesses can be attributed to the method of measurement 
required. 

Figure 3-2 shows the side view of anode plate #1 which is 
the view seen when looking at the end of the battery. All 
plates have essentially the same dimensions. The area of the 
plates is approximately 29cm 2 . 
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COMPONENT DISTANCE THICKNESS 

(mm) (mm) 






Anode Plate #1 


1.5 


1 .5 


AVERAGE 


THICKNESS 


Separator #1 


4.0 


2.5 


Anode 


1 . 57 5mm 


Cathode Plate #1 


6.2 


2.2 


S. Dev. 


0 .222 


Separator #2 


9.0 


2.8 


Var . 


0.0369 


Anode Plate #2 


10 . 3 


1.3 






Separator #3 


13.0 


2.7 


Cathode 


2 . 07mm 


Cathode Plate #2 


15.0 


2.0 


S. Dev. 


0.115 


« 

Separator #4 


17.7 


2.7 


Var. 


0 .0089 


Anode Plate #3 


19.5 


1.8 






Separator #5 


22.0 


2.5 


Separators 


2 . 6 2mm 


Cathode Plate #3 


24.0 


2.0 


S. Dev. 


0.133 


Separator #6 


26.5 


2.5 


Var . 


0.015 


Anode Plate #4 


28.2 


1.7 







Table 3-2. Measurements of the Distances of the Components 
from Edge of Anode Plate #1 and the Thicknesses 
of the Anode and Cathode Plates and Separators . 
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Figure 3-3. 



hole hole 




Table 3-3 shows the volume of the materials in a cell, 
including the deduced volume of electrolyte of 80cc. 

Access to the connectors between cells #1 and #2, and 
cells #3 and #4 was made by drilling into the top of the case. 
The location of the holes is shown in figure 3-3. 

The following observations were made of the battery in its 
unused condition. The anode plate viewed from the side had a 
light grey appearance; viewed from the end it was darker grey. 
Bubbles were seen on the plate and periodically several would 
rise to the surface. Cathode plates could only be observed 
from their ends due to the configuration of the cells. They had 
a dark brown rusty coating on their surfaces. Bubbles were not 
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INTERNAL VOLUME (per cell) : 

L = 29mm 
W = 51.5mm 

Height to upper electrolyte level (marked on case) = 80mm 
Volume = 120cc 

VOLUME OF ELECTROLYTE: 

Vol. of anode (per cell) = 29cm * .1575cm * 4 = 18.3cc 
Vol. of cathode (per cell) = 29cm * .207cm * 3 = 18.0cc 

Internal volume - plate volume = 120cc - 36cc = 84cc 
Assume solid volume of separators = 4cc 

VOLUME OF ELECTROLYTE = 80cc/cell 



Table 3-3. Volume of Materials within a Cell. 
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generally observed on these plates. 

These observations were made several weeks after the 
battery had first been activated, but before it had been 
used. When the battery had first been activated, the anode 
and cathode plates had more nearly the same appearance, which 
was light grey. 

The separators appeared to be made of a "fiber wool" 
placed on a paper (or felt) backing. The backing was always 
placed against the anode plates, while the "wool" was against 
the cathode plates. Small bubbles were entrapped in the fiber. 

The cavity below the plates had a very slight amount of 
rusty-colored, fine-grained sediment. 

3 . 2 EXPERIMENTAL SET-UP 

Figure 3-4 is a circuit diagram of the experimental set-up 
using a double pole-double throw switch. The circuit was 
configured in this manner so that no disconnecting and 
reconnecting would be required between charges and discharges 
and a single ammeter could be used in either mode. 

Two different battery charging set-ups were used. One was 
simply the use of a commercially available 6 Volt Tapering-Type 
battery charger. The battery charger was manufactured by 
Mallory, Indianapolis, Indiana (Engineering Projects Lab 
Number 1167), Model 6SAC10. 



Two chart recorders were used. One was a 4 channel chart 
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Figure 3-4. Circuit Diagram of the Experimental Set-Up 
Using a Double Pole-Double Throw Switch. 

recorder manufactured by Gould Brush 2400 Model 2007-4404-00 
Serial 00473. All four channels were used to record the 
voltage on each cell and the current simultaneously. This 
provided the experimenter with a complete picture of the 
behavior of the battery. 

The single channel chart recorder, a Hewlett-Packard 
Mosley 680 Strip Chart Recorder, Model H01-680, Serial # 
712-01825 (Engineering Projects Lab #1175) was used to record 
the voltage of the cell with the lowest capacity, i.e., the 
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cell which controlled the end of discharge point. This was 
cell #3 in battery #2. The single channel chart recorder was 
used twenty-four hours a day during life cycle testing, 
primarily because its slowest speed of 1 inch per hour gave a 
concise record of charge/discharge cycling. 

The four channel chart recorder was used when a more 
complete picture of the charge/discharge cycles was desired. 

The slowest speed on this chart recorder was 5mm per 100 sec 
(approximately 7 inches per hour) . Continuous running would 
have produced an unmanageable amount of paper. 

The set-up for life cycle testing included a sensing 
circuit to determine the end of discharge and end of charge and 
a relay to perform the switching of the battery circuit, as 
shown in figure 3-4. Figure 3-5 shows the schematic of the 
sensing and relay-driving circuit. The sensing circuit is a 
Schmitt Trigger described by Millman and Taub (1965) . 

It is a bistable circuit, and for a loop gain greater 
than 1, its input-versus-output displays hysteresis, as 
shown in figure 3-6. 

The relay driving circuit simply consists of a single 
transistor connected in series with the relay, with the out- 
put of the Schmitt Trigger to the base of transistor Q3 . The 
diode connected across the relay was simply for the purpose of 
providing positive opening of the relay when the driving 




Figure 3-5. Schematic Diagram of the Schmitt Trigger 

t 

Sensing and Relay-Driving Circuit. 

transistor turns off. The other relay shown in the circuit was 
included as a fail-safe device. It is energized by 115V ac 
line power. Its terminals are connected in series with the 
ammeter so that if there is a loss of power, the battery cir- 
cuit is opened. The terminals of the main relay were wired to 
the terminals of the double pole-double throw switch through a 
plug. When it was desired to perform controlled testing of 
the battery, the sensing and relay circuit was unplugged. 

One recommended modification to the circuit is to insert a 



3 position switch to the base of the driving transistor Q3 . 
One position would connect the Schmitt Trigger to the tran- 
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Figure 3-6. Hysteresis of Schmitt Trigger Circuit. 

sistor, as shown in figure 3-5. Another would bias Q3 
off, placing the battery on discharge, while the third pos- 
ition would turn Q3 on placing the battery on charge. This 
would facilitate making the controlled charges and discharges. 
Also, a simple open or closed switch in the fail-safe relay 
line might be added. With the set-up as shown , it was 
necessary to unplug the relay's energizing power from the wall 
receptacle . 

The Schmitt Trigger was very sensitive to fluctuations of 
short time duration. This might be improved by inserting small 
capacitors from the base of Q1 and Q2 to their emitters. 
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3.3 TESTING PROCEDURE 

There were basically two objectives of the testing. 

One was to experimentally obtain a family of discharge curves 
for the verification or modification of the Shepherd (1965) 
equation. The other was life-cycle testing in order to 
determine the change in battery characteristics with age. 

To obtain a family of discharge curves, constant current 
discharges were performed at 0.6, 1.5, 3.6, and 5.4A. These 
values were chosen because it was expected that they would 
provide a range of values similar to those shown in a plot 
of typical curves, such as those in figure 1-1. Once two runs 
have been made at different currents, Peukert ' s (1897) 
equation can be used to determine the current required for the 
discharge times desired. 

In an effort to place the battery in identical conditions 
prior to a discharge in which the data would be used for 
comparisons, the following procedure was followed. The 
battery was fully charged, followed by a pre-conditioning 
discharge at 3.6A. The reason for the pre-conditioning dis- 
charge was to remove variations in the "memory effect" (dis- 
cussed in Section 1.5) caused by previous discharges. The 
pre-conditioning discharge was followed by a constant current 
charge of 3.0A. The constant current was maintained by the use 
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Figure 3-7. Typical Charging Curve at Constant Current 
of 3.0A. 

of a 12 Volt battery connected in series with a variable 
resistor and was used as a charger. 

These values were chosen because they were within the 
range of currents used during life cycle testing and for con- 
venience. By using this scheme, it was felt that repeatable 
results could be obtained. The end of charge was determined to 
occur when the battery was gassing heavily and the cell voltage 
had remained constant for about 10 to 15 minutes, as shown in 
figure 3-7. 

Prior to discharging the battery, it was left on open 
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circuit for a period of at least one hour but not generally 
more than twelve hours. This would allow the electrolyte to 
become fairly uniformly distributed throughout the cell and 
also for the temperature to stabilize. 

During discharge, the variable resistor which served as 
the load was adjusted periodically to maintain constant 
current. The end of discharge occurred when the voltage of a 
cell dropped off abruptly. Discharge was terminated when this 
voltage dropped to approximately 1.0V. 

Life cycle testing was performed using the Schmitt 
Trigger circuit and relay with the Mallory battery charger 
and a constant load resistance of about 1.6ft. The charging 
current was in the vicinity of 3.0A, ranging from 3.5 at the 
start of charge to about 2A at the end. The average discharge 
current was in the vicinity of 3.5A. The Schmitt Trigger 
circuit was set to end discharge and commence charging when the 
battery voltage dropped to about 5V, and to switch on the high 
end at about 7.5V. The setting at the low end could vary by as 
much as several tenths of a volt without any particular change 
in the time of switching, due to the very rapid fall off of 
voltage at this point. The setting at the high end needed to 
be much more precise due to the leveling off of voltage at the 
fully charged point. This proved to be most difficult and 
introduced variations in the state of charge at the commence- 
ment of discharge under life cycle testing. 
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CHAPTER 4 

RESULTS AND DISCUSSION 
4.1 MODIFICATIONS TO THE SHEPHERD EQUATION 

The numeric data corresponding to figures 4-1 through 4-36 
in this section is found in appendix A. Figure 4-1 shows the 
family of discharge curves obtained following the methods out- 
lined in Section 3.3. Figure 4-2 shows the same curves after 
manual smoothing. One of the primary effects of this smoothing 
was the removal of the initial rapid drop of voltage, which 
Shepherd accounted for in his activation polarization term, and 
the "coup de fouet." Since the initial drop was completed 
within seconds of starting the discharge, and the "voltage dip" 
had disappeared within about two minutes, ignoring these pheno- 
mena should have little effect on the applicability of a model. 

The first step toward the establishment of a model was an 
attempt to fit the discharge data to Shepherd's equation: 




The values found by a least squares fit were: 

E = 2.295 
s 

K = 0.08086 
Q = 6.844 
R 0 = 0.00092 

with a correlation Z (E - E ) 2 for all curves = 3.5008. 

0Xp 
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4-1. Family of Experimentally obtained Discharge 
Curves. 
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Figure 4-2. Smoothed Family of Discharge Curves. 
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Figure 4-3 shows the resulting family of discharge 
curves. Comparisons of each individual discharge curve and 
its corresponding modeled discharge curve are shown in figures 
4-4 through 4-7. 

These results illustrate the need for some modification. 
The rapid drop-off of voltage at the end of discharge was only 
predicted for the .6 amp discharge. This was the expected 
result, considering the denominator of the concentration 
polarization term. 



2.4 
2.2 
2.0 
1.8 

0 1.6 

1 1-4 
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012 3 4567 
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Figure 4-3. Family of Curves for Curve Fit to Shepherd's 
Equation . 
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ampesie kouA-i 

e 4-4. 0.6 Amp Discharge Curve and the Modeled Curve for 

the Fit of the Family of Curves to Shepherd's Equation. 




ampe/ie houfU> 



Figure 4-5. 1.5 Amp Discharge Curve and the Modeled Curve for 

the Fit of the Family of Curves to Shepherd's Equation. 
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Figure 4-6. 3.6 Amp Discharge Curve and the Modeled Curve for 

The Fit of the Family of Curves to Shepherd's Equation. 
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Figure 4-7. 5.4 Amp Discharge Curve and the Modeled Curve for 

the Fit of the Family of Curves to Shepherd's Equation. 
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Next, a curve fit for each individual curve was per- 
formed using the Shepherd equation. Table 4-1 lists these 
coefficients with the graphical results shown in figures 4-8 
through 4-11. 



Current, i 


0.6 


1.5 


3.6 


5.4 


E 

S 


2.18 


2.18 


2.18 


2.178 


K 


0.0137 


0.0238 


0.0031 


0.00142 


Q 


6.502 


5.302 


4.373 • 


3.991 


Ro 


0.2779 


0.0481 


0.0734 


0.05924 


Correlation 


0.0718 


0.1239 


0.1585 


0.2638 



Table 4-1. Curve Fit of Shepherd's Equation to Individual 



Discharge Curves. 




ampete kou.A.4 

Figure 4-8. 0.6 Amp Discharge Curve and the Modeled Curve for 

a Fit to the Shepherd Equation of this curve alone. 
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re 4-9. 



1.5 Amp Discharge Curve and the Modeled Curve for 



a Fit to the Shepherd Equation of this Curve Alone. 
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Figure 4-10. 3.6 Amp Discharge Curve and the Modeled Curve for 
a Fit to the Shepherd Equation of this Curve Alone. 
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Figure 4-11. 5.4 Amp Discharge Curve and the Modeled Curve for 

a Fit to the Shepherd Equation of this Curve Alone. 
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If the slope of the model were adjusted then the cor- 
relation between discharge and model would be very good. This 
suggests modifying the resistance polarization term such that 
Vo = V 0 (q,i) = (aq+b)i. R 0 = aq+b takes into account that the 
ohmic resistance of the cell changes as the discharge reaction 
proceeds toward completion. This model simply predicts a 
linear change in resistance. It may be found that R 0 is also 
a function of current since the PbSOi, build-up occurs in a 
different manner at different discharge rates. 

The value for Q follows the Peukert equation very well 
when C and N have the values 5.803 and 1.2227 respectively. 
Figure 4-12 shows a log-log plot of Q versus I for Peukert' s 
equation . 
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Figure 4-12. Log-Log Plot of Peukert' s Equation 

( 1 -n) 



Q = Cl 



c = 5.803 n = 1.2227 
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Modifying the Shepherd Equation by the inclusion of 
Peukert's capacity model yielded a fit which predicted very 
well the drop in voltage at the end of discharge. The 
results obtained did not produce a very good fit to the left 
of this point, however. 

These plots are shown in figures 4-13 through 4-17. 
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Figure 4-13. Family of Curves for Modification of Capacity 




Term. 
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4-14. 0.6 Amp Fit with Capacity Term Modified. 
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Figure 4-15. 1.5 Amp Fit with Capacity Term Modified. 
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4-16. 3.6 Amp Fit with Capacity Term Modified. 
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Figure 4-17. 



5.4 Amp Fit with Capacity Term Modified. 



As mentioned in Section 2.2, the product of K and i, 
found from the curve fits figures 4-8 through 4-11, tended to 



be constant, as shown in table 4-2. 



i 


0.6 


1.5 


3.6 


5.4 


K 


0.0137 


0.0238 


0.00310 


0.00142 


K*i 


0.00822 


0.0357 


0.0112 


0.00767 



Table 4-2. Comparison of K and the Product K*i. 

This suggests that in equation 2.5, ^ V d (q,i), but rather 

= K q S q • B Y including both this new and 

capacity terms, the plots of figures 4-18 through 4-22 were 
produced. These show a definite improvement over previous fits 
for the family of discharge curves. The values determined were 



E = 2.002 
s 

K = 0.009 
R 0 = 0.03 



with a correlation of 1.8309. 



62 



2.2 

2.0 

1.0 

• 8 1.6 

<c 

5 14 

o ; • 

> 

1.2 

1.0 
0.0 

■0 1234567 

ampeAe houJiA 




Figure 4-18a. Experimentally Determined Family of Curves 

Modification of = V^(c) and Capacity are shown below: 
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Figure 4-18b. Family of Curves for Modification of 
= V^(q) and Capacity. 
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e 4-19. 0.6 Amp Fit with Modification of V^ = V^ (q) 

and Capacity. 
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Figure 4-20. 1.5 Amp Fit with Modification of V^ = V^(q) 

and Capacity. 
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Figure 4-21. 3.6 Amp Fit with Modification of = V^(q) 

and Capacity. 
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Figure 4-22. 5.4 Amp Fit with Modification of = V^(q) 

and Capacity. 
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Each plot of figures 4-18 through 4-22 suggests that if 
the ohmic resistance were modeled as a linear function of 
accumulated charge, a much better fit could be obtained. 

Figures 4-23 through 4-26 show the results of individual 
fits for each discharge curve by modeling 

Ro = Ro (q) = R a q + R^ . 

Table 4-3 lists these results. 



i 


0.6 


1.5 


3.6 


5.4 


E 

S 


2.180 


2.180 


2.180 


2.180 


K 


0.00744 


0.03453 


0.00876 


0.00615 


R a 


0.1183 


0.00073 


0.01881 


0.01771 


R b 


0.03867 


0.05424 


0.03253 


0.02016 


Correlation 


0.05279 


0.12685 


0.03129 


0.08875 


Figure 


4-23 


4-24 


4-25 


4-26 



Table 4-3. Curve Fit for Individual Discharge Curves with 
Modifications to Q, V^, and R 0 . 



VOLTAGE 
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Figure 4-23. Fit of 0.6 Amp Alone for Modification of 
Ro and Capacity. 
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Figure 4-24. Fit of 1.5 Amp Alone for Modification of 
Ro and Capacity. 
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Figure 4-25. Fit of 3.6 Amp Alone for Modification of 
R 0 and Capacity. 
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Figure 4-26. Fit of 5.4 Amp Alone for Modification of 
Ro and Capacity. 
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Tw© more curve fits to the family of discharge curves were 
performed using the modified ohmic resistance term and 
Peukert ' s capacity model. 

The first was performed using Shepherd's term with the 
results shown in figures 4-27 through 4-31. The second fit 
was made using the modified term. These results are shown 
in figures 4-32 through 4-36. 
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Figure 4-27a. Experimentally Determined Family of Curves 

Modification of Capacity and R 0 = Ro (q) are shown 
below: 




ampzfiz houA-i 

Figure 4-27b. Family of Curves with Modification of Capacity 
and R 0 = R 0 (q) • 
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4-28 0.6 Amp Fit for Modification of R 0 and Capacity. 




Figure 4-29. 



1.5 Amp Fit for Modification of R 0 and Capacity. 



VOLTAGE £ VOLTAGE 



77 




01234567 

ampere kou-U 



4-30. 3.6 Amp Fit for Modification of Ro and Capacity. 




Figure 4-31. 5.4 Amp Fit for Modification of Ro and Capacity. 
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Figure 4-3 2a. Experimentally Determined Family of Curves 

Modifications to V^, Ro , and Capacity are shown below: 
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Figure 4-32b. Family of Curves with Modifications to 
V d , Ro, and Capacity. 
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e 4-3 3. 0.6 Amp Curve with Modifications - to V d , R 0 , 

and Capacity. 
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Figure 4-34. 1.5 Amp Curve with Modifications to V^, R 0 , 



and Capacity. 



VOLTAGE 




ampere /ioua<6 



Figure 4-35. 3.6 Amp Curve with Modifications to V d , R 0 

and Capacity. 
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Figure 4-36. 5.4 Amp Curve with Modifications to V^, Ro 

and Capacity. 
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4.2 LIFE CYCLE TESTING 

As was previously mentioned in Section 1.5, according to 
Bode (1977), the ideal requirement for life testing is that the 
life test correspond to the load conditions of the battery in 
the field. However, since ideal testing is impossible due to 
conditions impossible to duplicate in the laboratory, one 
must choose another type of test. Kordesch (1977) states that 
the Cycling Life Test is a test principally to evaluate battery 
internal components. It does not attempt to simulate car 
service conditions which are extremely variable. Naturally, 
this is the type of test chosen for testing the lead-acid 
battery which was used for experimentation for this thesis. 

The life-cycle testing procedure used to determine the 
change in battery characteristics with age is outlined in 
Section 3.3. Some of the pertinent data concerning life- 
cycle testing follows: 



Number of Cycles = 135 
Time on Charge = 193 hours 
Time on Discharge = 114 hours 

Total number of Ampere Hours Discharged = 375 
Amount of Water Lost Due to Hydrolysis = 12ml/cell 



To determine changes in cell capacity during the life 
of the battery several controlled discharges at 3.6 amps were 
made at different times during the testing. Table 4-4 
shows these results. The final capacity was 94% of the pre- 
vious capacities. Further testing is required to determine 
the rate at which the capacity decreases . 



Number of 


Prior Cycles 


25 


91 


135 


Discharge 


Time 


72 


72 


67.5 


Amphere-hour Capacity 


4.32 


4.32 


4.05 



Table 4-4. Capacity Change of Cell #3 During the Life 
of the Battery. 

A great deal of data was generated during the testing but 
it was considered far to voluminous to be included in its raw 
form in this report. This data, which is primarily in the form 
of strip chart recordings, was placed on file with Professor 
H.M. Paynter at the Massachusetts Institute of Technology. 

4.3 UNCERTAINTIES IN EXPERIMENTAL RESULTS 

Uncertainties involved in the experimental results were 
due partly to imperfect knowledge of the state of charge and 
the inability, at least with the experimental set-up used, to 



control certain parameters. The most notable of these was 
temperature. Variations in temperature have been recognized 
to produce significant changes in battery characteristics . 

Great variations in temperature were not encountered since the 
experiment was performed at room temperature, but there were 
some variations in this temperature by nearly as much as 10°C. 
The effect of internal heating at higher discharge rates may 
have affected results. 

The difficulty of knowing the precise state of charge also 
added a degree of uncertainty. Included in this was the 
"memory effect," since the amount of preconditioning required 
to eliminate it was not known. 

Also, the knowledge of the mechanism of lead sulfate 
build-up at the various discharge rates was not available. 

Due to the differences of the results during curve fitting 
for the 1.5 amp discharge with respect to the other curves, 
there was some suspicion that the conditions for that particu- 
lar discharge may have been different than for the others. 

4.4 ADDITIONAL OBSERVATIONS 

The method of charging the battery and the point at which 
charging was terminated was deemed to be the most important 
factor in placing the battery in the same state-of-charge 
before the controlled discharges. Figure 4-37 shows a typical 
constant current charging curve for this experiment. 
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Figure 4-37. Typical Constant Current Charge. 

The amount of gassing in the battery was determined by 
visually observing the amount of gas evolution. The clear 
battery case was vital in making this determination. By 
observing the top of the plates from the side of the battery, 
the following amounts of gassing were determined. Slight 
gassing was determined when occasional bubbles rose to the 
surface. This was about the same as when the battery was 
standing. Light gassing occurred when there was a regular 
stream of these bubbles rising to the surface. Moderate gas- 
sing existed when considerably more activity was observed. 
This was accompanied on the charging curve by a noticeable 
increase in the slope. Heavy gassing was typified by a sheet 
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of small bubbles rising and a generally more violent condition 
at the surface of the electrolyte. This occurred when the 
slope of the charge curve reached its maximum and continued 
throughout overcharging when the voltage stopped increasing. 

To obtain equivalent charged conditions, the battery was 
overcharged when essentially no further increase in voltage 
occurred at constant current for about 15 to 30 minutes. Most 
of the energy into the battery at this point went into the 
hydrolysis of water, indicating that the charge reaction was 
essentially completed. The observation of the "coup de fouet" 
after this type of charge and its absense after a partial 
charge, adds some validity to the degree of completion of the 
reaction . 

An explanation of this phenomenon is that upon initiation 
of discharge, the electrolyte goes into a meta-stable state, 
super-saturated with PbSC >4 , until PbS 04 crystallization sites 
start to form on the plates. Once this occurs, the voltage 
recovers from its "dip" and proceeds in a smooth fashion. This 
recovery was observed to occur after about two minutes . The 
"coup de fouet" was not observed on a discharge following a 
partial charge, and was significant for discharges at higher 
discharge rates, 3.6 and 5.4 amps, which followed a full 
charge. The magnitude of this dip was approximately 0.01 volt. 

When the battery was new, Cell #2 had the lowest capacity. 
After several days of testing, while moving the experimental 
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set-up. Cell #3 was inadvertantly placed across a . IQ resistor. 
After that occurred. Cell #3 had the lowest capacity and Cell 
#1 had the greatest. 

After all testing was completed, the following observa- 
tions of the battery were made. The positive plates and 
connectors had a dark brown color (artist's color "Burnt 
Umber"), which had a fairly rough appearance. There was some 
spalling on one edge of one of the plates in cell #3. 

The negative plates had a battle-ship grey color. One 
small piece of material was observed to have broken away from 
the plate. 

There was a fairly uniform coating of sediment in the 
bottom of the battery. This sludge consisted of two distinct 
materials. One was a light-colored material with a slightly 
pinkish tint. The other which uniformly covered the bottom 
of the battery, and was layered on top of the other sediment, 
was a dark brown color much the same as the color on the 
positive plates. The light-colored material may have been 
white rather than slightly pink because it had been observed 
by the author et al that the battery case itself seemed to cast 
the pink tint onto the interior components. It is believed 
that the light-colored material is actually predominantly 
PbSO 4 . 

Vinal (1955) states that overcharging produces sludging of 
a fine brown sediment, and that prolonged undercharging causes 
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a fine white powder which is principally PbSOt, to be deposited 
on the bottom. It is felt that the dark brown sediment is 
Pb02 1 which represents a loss of active material from the 
positive plates. The author notes that prolonged undercharging 
took place during life-cycle testing. 



CHAPTER 5 



CONCLUSIONS AND RECOMMENDATIONS 



The Shepherd (1965) equation 

E = E s " V d (i,q) - Vo (i) - V a (q) (5.1) 

V (q) = activation polarization = A exp ^2 (5.2) 

a v 

(neglected due to the short time duration) 

Vo (i) = ohmic polarization = Roi (5.3) 

V d (i,q) = diffusion polarization = K ^ ^ — i (5.4) 

proved to be an excellent starting point for establishing an 
electrochemical state model for lead-acid batteries. Using the 
lead-acid battery of this investigation, the author was able to 
get a fit to individual discharge curves which required only 
the modification of the ohmic polarization term (equation 5.3) 
to the form 

V 0 (i,q) = (R a q + R b ) i. (5.5) 



in order to get an excellent fit. The results obtained from 
this fit suggested that a better fit for the family of curves 
might be obtained by setting K*i constant. 



The resulting modified term was 




(5.6) 



This was the form of the transport polarization suggested by 
Bode (1977) . They also show that the use of the Peukert 
(1897) capacity model 



would predict the rapid drop of voltage at the end of discharge 
very well. 

The modification of the and Q terms seemed valid 

modeling for the family of discharges, but difficulties were 
apparent with the ohmic polarization term. There were no 
particular consistencies in the values of R & and for the 
different discharge rates. This may be due to the fact that 
the buildup of PbS0 4 at the plates occurs in a different manner 
for different discharge rates. Bode (1977) shows that for 
discharge rates of 6 mA/cm 2 and 1 mA/cm 2 the distribution of 
PbS0 4 throughout the positive, capacity limiting plate, is 
uniform at full discharge. For current densities of 30 mA/cm 2 
and 180 mA/cm 2 , there is a lower PbS0 4 content at the center of 
the plate while the content near the surface is the same as 



Q (I) = CI (l-n) 



(5.7) 
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for the lower current densities. 

The cell under test had current densities of 3.4, 8.6, 
20.7, and 31 mA/cm 2 for discharges of 0.6, 1.5, 3.6, and 5 . 4A 
respectively. This suggests that the same mechanisms were 
encountered in the battery discharge. Since this effect is 
clearly a function of current, a modification of V 0 (i,q) other 
than that suggested in equation 5.5 may yield a better fit. It 
may have a form similar to Shepherd's original term, 
equation 5.4. 

Another consideration might be that since equation 5.4 was 
based on a small current density approximation, another look 
might be made toward applying equations 2 . 7 and 2 . 8 upon which 
the approximation was based. 

A very exciting area of future endeavor, and one which 
could prove very rewarding, would be the establishment of a 
model which would predict the family of discharge curves. 

Tremendous advancements in the lead-acid battery have been 
made since the time of Plante and this progress should continue 
steadily due to the importance of this device in myriad diverse 
applications. From the standpoint of the System Designer, the 
establishment of a good state model is a very exciting pros- 
pect. This would then allow for a determination of the optimum 
use of the battery, subject to the constraints of the applica- 
tion. 

Further effort in the direction of establishing a model 



which would predict the family of discharge curves should 
include the effect of discharges at various rates during a 
single cycle. Other effects such as pulsed discharge including 
both the frequency and duty cycle could then be added. 

The inclusion of the prediction of the change in internal 
characteristics during the secular life of the battery is 
virtually essential for an optimization. The objective in the 
design of almost any system is the optimization of a cost 
(possibly acquisition cost or life cycle cost) subject to 
various operational criteria. For a system such as an electric 
vehicle, a very significant cost is the repair cost of battery 
replacement. Therefore, it is necessary to be able to predict 
the lowest cost scheme for charge and discharge battery requir- 
ments and replacement frequency for various charge/ discharge 
schemes. To this end, further life cycle testing is recommend- 
ed. 

To produce applicable results, it is felt by the author 
that all cycles should be identical except those due to secular 
changes in the battery. A rather sophisticated control system 
needs to be developed. All charges and discharges should be 
conducted at constant current. During charge, a significant 
period on overcharge is recommended to ensure that the state- 
of-charge is nearly identical prior to any discharge. Since 
the amount of overcharge cannot be determined by sensing the 
voltage level alone, some means needs to be devised to deter- 
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mine the slope of the voltage curve and to control the time 
on charge at this level- It is also recommended that a single 
cell be used during testing to eliminate any effects from the 
other cells. 

An attempt has been made throughout this thesis to draw 
conclusions on the suitability of the Shepherd (1965) equation 
as a state model for lead-acid batteries. Modifications to 
this equation were proposed, among them was the modeling of 
the ohmic polarization as a linear function of accumulated 
charge to obtain a good fit to a single discharge curve. These 
curves revealed the interesting result that the product K*I 
was nearly constant. Also, Peukert's equation was used to 
model the cell capacity. 

The difficulties encountered in the attempt to obtain an 
optimum state model have been evaluated and enumerated. How- 
ever/ new discoveries about the lead-acid battery are being 
made, and new technology may supply the means to overcome the 
difficulties. Almost certainly, economic, ecological, socio- 
logical and technological problems with our current energy 
resources will dictate to society the necessity of assessing 
the place of the lead-acid battery in our future. 
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APPENDIX A 



NUMERIC DATA CORRESPONDING TO FIGURE k-1. 



0 . 6 Amp 


Di s ch 


Ah 


Volt 




( exp ) 


0 • 0 0 


2.110 


1*80 


2.060 


3*80 


2.010 


a . 2f. 


1 .<590 


*.80 


1 .970 


'5.4 0 


1 .910 


5 * 7 C 


1 .830 


8.0 0 


1.830 


8*12 


1.810 


8.2* 


1 . 770 


: 8 * 3 0 


1.740 


8.38 


1.690 


.8 * 3 0 


1.63 0 


'8.4 2 


1 .470 


6.44 


1.030 



1.5 Amp Disch 





Volt 




( exp ) 


r . 00 


2 . 0 S 0 


C .75 


2 .- 0^0 


1.50 


2.020 


2.25 


2.000 


2*83 


1.990 


3 . 0 C 


1 . ° 3 C 


3.3 8 


1 .93 0 


4.13 


1.910 


4.50 


1 . 8 S 0 


4.85 


1 . 8*0 


.4.80 


1.610 


4.88 


1.770 


4.55 


1.720 


5.03 


1.600 


5.10 


1.100 


5.13 


0.85 0 



3.6 Amp Disch 



.h 


Volt 




( exp ) 


r .00 


2.00 0 


0.30 


2 • 0 1 0 


0.8 0 


2.00C 


0.9 0 


1.990 


1.20 


1.970 


1.50 


1 . c 8 0 


i . 8 0 


1 .94 0 


2.10 


1.930 


2.4 0 


1.910 


^ *r n 

C • 1 1/ 


1 .990 


3*00 


1 . 3 8 n 


3 * 3 0 


1.830 


3 * 0 0 


1 . 79 f' 




1 . 75 C 


5*98 


1.710 


4.08 


1.850 


4 . 1 n 


1.590 


4.20 


i . 4 0 1 


4.28* 


1 . 2 8 0 


4.3 2 


'O . 94 Q 



5 . k Amp Disch 
Ah Volt 

( exp ) 



— — — 


- 


— — 


0.00 


1 


• ^9 0 


9 H . 


p 


. 0 0 0 


o . 9 :• 


x 


. c 7 0 


1.35 




• 940 


1 * S 0 


1 


. 92 0 


2.25 


1 


.8 9 0 


2 . 7 0 


1 


. * 5 " 


3.15 


1 


♦ 79 0 


3.4 2 


1 


*79 3 


T . <*•? 


1 


.700 


!" . c 


1 


• 88 1 


3.78 


1 


. 5 7 0 


3.87 


1 


.230 


7 .96 


U 


. 0 0 C 
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Numeric data corresponding to Figures k-k and U-5 



Coefficients for the Model 

E s K Q R 0 Corr 

2.295 0.08086 6.84U 0.0092 3.5008 



0.6 Amp Discharge 




Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 




o.nc 


2.110 


2.246 


1 ^ 
1 ^ 
1 ^ 
| . 
1 ° 


1.80 


2.060 


2 • 2 2 R 


G. If 8 


3 . 6 C 


2.010 


2.192 


0.182 


4.20 


1 . n 9 9 


2.169 


0. 1 79 


4.80 


1.870 


2, . 1 0 2 


C . 1 6 2 


5.4 0 


1.810 


2.064 


0.154 


5.7 0 


1 . 8 8 G 


3.004 


0.1 24 


6.0 0 


1.330 


1 . 901 


C . n 7 1. 


6. 1? 


1.810 


1.836 


0.026 


6.2 4 


1 . 77D 


1.745 


-0. 025 


6.30 


1.740 


1.684 


-0.056 


6.36 


1.680 


1.608 


-0. 182 


6.38 


1 .630 


1 . c 63 


-0.067 


6.42 


1.470 


1.511 


0.041 


6.44 


1 . 030 


1.473 


0.443 



1.5 Amp Discharge 



Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 




1 O 
1 

1 • 
1 


2.059 


? • 1 72 


G . 1 20 


0.75 


2.040 


2.157 


C . 1 1 7 


1.5 0 


2.02 C 


2.138 


0.118 


2.25 


2 * 0 0 0 


2.113 


C . 1 1 3 


2.63 


1.990 


-.007 


0 . ] (17 


3. 0 0 


1 .98C 


2.077 


0 . : J 7 


3*38 


1 . 963 


2 • 09 4 


0 . ! : ' 1 


4.13 


1.920 


1 . “• .? 8 


0 . C 68 


4.^0 


1.880 


1.939 


0 . :• r 9 


4.6 5 


1.85 0 


1 .-15 


r .063 


4 • 8 C 


1.610 


1.887 


C . C 77 


4.98 


1 . 77 0 


1 .872 


r « *1 0 

- • i l c. 


4.95 


1.720 


1.653 


0.1; 5 


5 . C 3 


1.600 


1 . S3 7 


0 . 2 7 7 


5.10 


1.100 


1.817 


C . 7 17 


5.13 


0.850 


1 . 8 1 0 


C . 0 t- 0 
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Numeric data corresponding to Figures 4-6 and 4-7. 



Coefficients for the Model 

E s K Q R Q Corr 

2.295 0.08086 6.844 0.0092 3.5008 



3.6 Amp Discharge 




5.4 


Amp Discharge 




Ah 


Volt 


Volt 


Diff 


Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 






(Exp) 


(Eqn) 




0.00 


2.017 


2.000 


-0*017 


0 . C G 


2.014 


1.853 


-0.2 61 


G • 3 0 


2.010 


1 .987 


-0*023 


0.45 


1.590 


1 .92? 


- r . : c p 


0 .6 G 


2.000 


1 . Q 72 


-0* 02* 


0.90 


1.971 


1 « 7 87 


- C ♦ 1 h 3 


0.90 


1 .984 


1.956 


-G . 0 28 


1 .35 


1 .95 0 


1.746 


-0.204 


1.20 


1 .970 


1.938 


-0 . 0 32 


1.80 


1.920 


1.697 


- 0 . 2 r 3 


1.50 


1.960 


1.919 


-0*041 


2.2 5 


1.893 


1.939 


- n . 2 s 1 


1.8G 


1 .98 5 


1.896 


-0 * 0 49 


2.70 


1 .850 


1.569 


-0.261 


2.10 


1 .930 


1.871 


-0. G 59 


3.15 


1.800 


1 .481 


-0.31° 


2.40 


1.910 


1 .843 


-0* 0 67 


3.4 2 


1.750 


1.417 


-C . 3 ” 


2.70 


1.890 


1.811 


-0. 079 


3.60 


1.700 


1 .5 69 


- 0 . T '• 1 


3.0 0 


1 .86 0 


1.773 


-0.087 


3.69 


1.660 


1 .542 


- C - •: 2 8 


3.30 


1.830 


1.729 


-0.101 


3.78 


1.570 


i .314 


- 0 • 9r 6 


3 . 6 C 


1.790 


1 .677 


-0.113 


3.8 7 


1.230 


1 • 7 h 5 


C. 


3.78 


1.760 


1.641 


-0. 1 19 


3.93 


0 . fi G 0 


1.254 


G .454 


3.96 


1.710 


1.601 


-0.109 










4.08 


1.650 


1.571 


-0. Q 7? 










4.14 


1 .590 


1.555 


-0.035 










4.20 


1.480 


1.558 


0 . G5P 










4.26 


1.280 


1.520 


C .240 










4.32 


0.940 


1.502 


0.562 
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Numeric Data corresponding to Figures H-8 and U — 9 • 



Coefficients for the model. 



See table 4-1 



0.6 Amp Discharge 




Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 




0. 00 


2.110 


2.005 


-0.105 


i .ee 


2 • 0 fc 0 


2 .002 


- C . C 58 


3.6 0 


2.010 


1 . c 95 


-0.015 


4.2C 


1 . 990 


1 . 990 


0.0^0 


* . 8 0 


1.970 


1 .982 


0.012 


5.4 0 


1.910 


1 .965 


0 . 0 55 


5. 7C 


1 .9S0 


1 . <M7 


0. 0*7 


6.00 


1.830 


1 .907 


0.C 77 


6.12 


1.813 


1.873 


0. 0 63 


b.24 


1 .770 


1.809 


0 . C 39 


6.30 


1.740 


1 .748 


0.008 


6.36 


1.690 


1 .636 


-0 . C 54 


6.39 


1 .630 


1.535 


-0.0 55 


6.42 


1 .47 0 


1.361 


-0 . 1 09 


6.44 


1.030 


1.15 0 


C . 1 20 



1.5 Amp Discharge 



Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 




0 .CO 


2 .05 2 


? . '■ 7 *> 


o 


0 .75 


2.040 


2.060 


o . 


1 .50 


2.0 20 


2.0 5' 


C . C 3 ° 


2.25 


2.000 


'.04' 


o . : 1 5 


2.63 


1.990 


2 . 0 3 7 


: . 7 


3.0 0 


1.98 C 


2 . 02 * 


C. 4 5 


3.3° 


1.96* 


2 . 01 ' 


0 . 0 't 7 


4.13 


1.920 


1 . = 4 7 


0.027 


4.50 


1 . 8*0 


1.872 


- 0 . 00 ° 


4.65 


1.550 


1.8 1 = 


- 0 . 0*2 


4.80 


1.910 


1.731 


C 

r 

i 


4 . 8 a 


1.77' 


1.665 


-0 . 1 * 5 


4 .95 


1.720 


1 . c 7 0 


-0 . - 5 0 


5.0 3 


1.60 0 


1.425 


- 0 .175 


5.10 


1.10 0 


1.171 


0 . " 7 1 


5.13 


0 . S 5 0 


1.03 3 


0.1*9 
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Numeric data corresponding to figures U-10 and 4-11 



Coefficients for the Model 



See Table 4-1 



3.6 Amp Discharge 5.4 Amp Discharge 



Ah 


Volt 


Volt 


Diff 


Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 






(Exp) 


(Eqn) 




o . c c 


?.cr 


1.890 


-0.127 


0.0 0. 


2.014 


1 . 0 60 


-0 . 1 fc 4 


0 .30 


2.010 


1.389 


-0.121 


0.45 


' .9 9 0 


1. 9 «i 


-0.141 


0 o 


2.000 


1 . 3 P 3 


-0.112 


0.90 


1 ,97C 


1. s 4«. 


-0 .122 


0.90 


1.934 


1.887 


-0.09’’ 


1.35 


1.950 


1.84V 


-0. V 3 


1 .20 


1.970 


1.885 


- 0 . 0 ° 5 


1 .80 


1.920 


1 . 0 4 4 


-0.076 


1 .50 


1.950 


1.834 


-0.076 


2.26 


1.390 


1.04 1 


-0 . C 4 0 


1.30 


1.945 


1.93? 


-0 .('53 


2 .70 


1 . 8 c ■: 


1 . 8 3 4 


-o .r *6 


2.10 


1.9 30 


1.979 


-0.051 


3.15 


1 .EGG 


1.822 


0 .C->2 


2.4 0 


1.910 


1.576 


-0.034 


3 .42 


1.750 


1 . a 0 * 


C .' s 5 


2.70 


1.390 


1.872 


-0.019 


3 .60 


1 .70 0 


1.73 0 


0 . r c o 


3.0 0 


1.860 


1.365 


0.03 5 


3.69 


1.65' 


1.755 


0 . r 96 


3.30 


1.8 3C 


1.855 


0.025 


3 .7<> 


1.578 


1.713 


.14 3 


3.50 


1 .790 


1.338 


0.049 


3 .87 


1.23' 


' . c 0 5 


C . ?"5 


3.73 


1 .750 


1.819 


0 .059 


•7 .96 


0 . 8 r •" 


3 . 3 " 1 


C . 1 


3 .96 


1.710 


1.73 3 


0.073 










4.03 


1.550 


1.734 


3.034 










4.14 


1.590 


1.691 


0.101 










4.20 


1 . 480 


1 . .6 1 9 


0 .139 










4 .25 


1 .280 


1 .U5Q 


0.189 










4.32 


0.940 


0.980 


0 .040 
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Numeric Data corresponding to Figures t-lU and 4-15. 



COEFFICIENTS FOR THE MODEL 



E s 


K 


c 


n 


Co rr 




R o 




1.872 


0.00177 5- 


803 1.22270 


3 . 19808 




0 . 00651 


0.6 Amp Discharge 




1.5 


Amp Discharge 




Ah 


Volt 


Volt 


Diff 


Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 






(Exp) 


(Eqn) 




0 .00 


2.11° 


1.662 


-0.243 


0 . r, 0 


2.052 


1 . a 6 0 


-0.172 


1 .60 


2.060 


1.867 


-0.173 


0 .75 


2.04 r, 


1.°57 


-0.131 


3.60 


2.010 


1.366 


-0.144 


1 .50 


7.02 0 


1.050 


-0.161 


4.20 


1.990 


1.365 


-0.125 


2.25 


2 . 0 C * 


1 .°5 C 


-0 .142 


4 .80 


1.970 


1.364 


-0.106 


2 .6 3 


1.99*' 


1.°5' 7 


-0.133 


5.UC 


1.910 


1.662 


-0.048 


3 .0 0 


1 .98P 


1 . 8 5 5 


-0 . 1 24 


5.70 


1 .880 


1 .85 9 


-0.021 


3 . 3 a 


1.963 


1 . ? 5 5 


-0 . 1 A ? 


6 .00 


1.330 


1 .654 


0 . C 24 


4.13 


1 .92 A 


1.350 


-0 .070 


6.12 


1.310 


1.852 


0.040 


4 .50 


1 .3 80 


1 . a 45 


-0 .C a 5 


,6.24 


1.770 


1.342 


0 . C 7 2 


4.65 


1 .350 


1.341 


- 0 . C A 9 


'6.30 


1.740 


1.3 34 


0 .C a 4 


4 . c 7 


1 . 3 1 7 


1.834 


0 . A 7 4 


5 .36 


1.693 


1.812 


0.129 


4 .8 a 


1.770 


1 . n 27 


0 . A 5 9 


6 .30 


1.630 


1.80 6 


0.176 


4 .9 5 


1.720 


1.-22 


0. 102 


6 .42 


1.470 


1.734 


0.314 


5 . o' j 


1.600 


1.311 


0.711 


6 .44 


1.030 


1.757 


0.727 


5.10 


1 .107 


1.792 


0.692 










5.13 


0.850 


1.783 


0.933 
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Numeric Data corresponding to figures b-1 6 and 4-17 



COEFFICIENTS FOR THE MODEL 



E s 


K 


c 


n 


Corr 




R o 




1.872 


0.00177 5-803 1.22270 


3.19808 




0.00651 


3.6 Amp Discharge 




5.4 . 


Amp Discharge 




Ah 


Volt 


Volt 


Diff 


Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 






(Exp) 


(Eqn) 




0 .00 


2.01 7 


1.842 


-0 . 17 = 


0 .00 


2.014 


1.927 


-0.137 


0 .30 


2.010 


1 .342 


-0.163 


0 .4 5 


1 . 990 


1.326 


-0 ,1 c -4 


0.60 


2.000 


1.341 


-0.159 


C .90 


1.970 


1.325 


-0.145 


0.90 


1.984 


1 . 3. 4 1 


-0.143 


1 .35 


1.950 


1.922 


-0.12= 


1 .20 


1 .970 


1 . R 4 0 


-0.130 


1 .80 


1.923 


1.919 


-0.101 


1 .50 


1 .960 


1 . 0 3 9 


-0.121 


2.26 


1.894 


1.915 


- 0 . C 7 5 


1 .9 0 


1.945 


1.335 


-0 .107 


2 .70 


1.950 


1.907 


-0.043 


2.10 


1 .930 


1.335 


-0.044 


3.15 


1.300 


1.791 


-0.A9Q 


2 .40 


1.9 10 


1.334 


-0.076 


3 .42 


1.750 


1.770 


0 . G 2 9 


2.70 


1.890 


1.332 


-0.053 


3.60 


1.700 


1.73 s 


0.03 = 


3.00 


1.860 


1.828 


-0.032 


3.69 


1 . 6 6 r 


1.70 = 


0 . C 4 s 


3.30 


1.8 30 


1.822 


-G.09P 


3 .79 


1.57" 


1.652 


0.7 = 3 


3 .6 0 


1 .790 


1.812 


0.022 


3 .37 


1 .23" 


1 .50 = 


0 .2"9 


3 .73 


1.760 


1 .°01 


0.041 


3.95 


0 • 8 C 0 


0.972 


-0.429 


3.9 6 


1.710 


1.780 


0.070 










4.0 8 


1.650 


1.750 


0.100 










4.14 


1.590 


1.724 


0.134 










4.20 


1.480 


1.673 


0 . 198 










4.25 


1.28'' 


1 .57* 


0.298 










4 .32 


0 .940 


1 . 19"' 


0.257 
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Numeric Data corresponding to Figures 4-19 and 4 - 20 . 



COEFFICIENTS FOR THE MODEL. 



= s 


K 


c 


n 


Cor r 




R o 




> .002 


0.00900 5. 


803 1.22270 


1.83086 




0 .03006 


0.6 Amp Discharge 




1.5 Amp Discharge 




Ah 


Volt 


Volt 


Diff 


Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 






(Exp) 


(Eqn) 




0 .00 


2. Ilf' 


1.975 


-0.135 


0.00 


2.052 


1.948 


-0 . 1 '4 


1 .30 


2.050 


1 .972 


-0 . 0 0 0 


0.75 


? . 0 4 C 


1.947 


-0 .0 r *U 


3.60 


2.01" 


1.364 


-0 .045 


1 .59 


2.0 27- 


1 . 7 £• <4 


-0 .076 


4 .20 


1.990 


1 .959 


-0 .031 


2.25 


2.007 


1.341 


- C . C 3 9 


4.30 


1 .970 


1 .950 


-0 .020 


7.63 


1.957 


I.750 


-7 . 0 r 1 


5.40 


1.910 


1 .93 1 


0.021 


3 . C 0 


1.985 


1.336 


-0 .'44 


5 .70 


1 . S 8 0 


1.3 11 


0.031 


3.39 


1.963 


1.532 


-0.091 


6 .00 


1.830 


1 .867 


0 .037 


4.13 


1.920 


1.517 


-0 . 004 


6.12 


1.810 


1.931 


0 .021 


4.50 


1 . 8 p 0 


1.99"' 


0.017 


6.24 


1.77C 


1.351 


-0.009 


4.65 


1.930 


i.saii 


0.C34 


6.30 


1 .740 


1.694 


-0.046 


4 .80 


1.817 


1.867 


0.57? 


6.36 


1 . S 9 C 


1 .572 


-0.118 


4.33 


1 • 7 7 j 


1.84" 


0 .075 


6.34 


1 . 6 3 A 


1,451 


-0 .16 9 


4.95 


1.720 


1 . 5 2 1 


0.171 


6 .42 


1 .470 


1.270 


-0.200 


5 .03 


1 . 6 :■ r 


1 7QT 

• ■* 


0 . 1 5 5 


6.4 4 


1.030 


1.040 


0.010 


5.10 


1.10 0 


1.720 


0 . 6 00 










5.13 


n . 8 c C 


1 . ^ 3 7 


0.737 
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Numeric data corresponding to Figures 4-21 and 4-22. 



COEFFICIENTS FOR THE MODEL 



E s 


K 


c 


n 


Corr 




R 0 




2 .002 


0.00900 5. 


803 1.22270 


1.83086 




0 .03006 


3.6 Amp Discharge 




5.4 Amp Discharge 




Ah 


Volt 


Volt 


Diff 


Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 






(Exp) 


(Eqn) 




0.00 


2.017 


1.9 35 


-0.132 


0 .00 


2.014 


1 . 3 3 1 


-0 . 1 G 3 


0 .30 


2.01 0 


1 . ? 8 4 


-0 .126 


0.45 


1 .3 7 0 


1 .93 n 


-0 .1*0 


0 .60 


2.000 


1 . 9 87 


-0.117 


0.90 


1.970 


1 . 0 2 5 


-0 .142 


0.90 


1 .934 


1.382 


-0.102 


1 .35 


1.950 


1.92* 


-0.124 


1 .20 


1.970 


1.831 


-0.089 


1.80 


1.92C 


1.523 


-0.097 


1 .50 


1.960 


1.330 


-0 . Q80 


2 .25 


1.897 


1.810 


-0 .O^l 


1 .30 


1.945 


1.378 


-0 . 0*7 


2 .70 


1 .8 50 


1 .» 12 


-0 .038 


2.10 


1.930 


1 . 5 7 6 


-0.054 


3.15 


1 . 8 C n 


1.777 


-0.003 


2.4 0 


1.910 


1 . 0 74 


-0.036 


3.42 


1 . 7 5 r 


1.77* 


r\ a c. 

J • . t 


2.70 


1 .390 


1 . 3 70 


-0.020 


3 .60 


1 .70: 


1.747 


0 . 0 u 7 


3.00 


1.860 


1.355 


0.005 


3 .go 


1.55' 


1.719 


0 . r C P 


3.30 


1.830 


1 . «57 


0.027 


2 .78 


1 .5" 7 -' 


1.56* 


r. 

v • 


3.60 


1 .79 


1.842 


0 . 0 C 2 


3 .67 


1 .23: 


1.53C 


0 .370 


3 .73 


1.760 


1.526 


0.066 


3 .96 


0 . 5 0 0 


0.460 


-0.340 


3 .96 


1.710 


1.79 6 


0.0O6 










4.03 


1.650 


1.755 


0. 105 










4.14 


1.590 


1.717 


0.127 










4 .27 


1 .450 


1.55? 


0 . 172 










4.25 


1.230 


1.511 


0.291 










4 .3 2 


0.940 


0.974 


0.034 
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Numeric data corresponding to Figures 4-23 and 4-24. 
COEFFICIENTS FOR THE MODEL 
SEE TABLE 4-3. 



0.6 Amp Discharge 



Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 




0 .0 0 


2.110 


2.102 


-0.0 08 


1.30 


2.060 


2.G57 


-0.0 03 


3.6 0 


2.010 


2.009 


-n . o oi 


o 

CM 

• 


1.990 


1 .991 


0.0 01 


4.6 0 


1.970 


1 .969 


-0.001 


5.40 


1.910 


l. Q 40 


0.0 30 


5.70 


1.830 


1.916 


0. 0-36 


6.00 


1.830 


1.373 


0.043 


6.12 


1.810 


1.8*0 


0.0 30 


6.24 


1 . 770 


1 .780 


C . 0 1 0 


6.3 0 


1 .740 


1 .723 


-0.017 


6.36 


1.690 


1 .621 


-0.069 


6.39 


1.63G 


1 .529 


-0.101 


6.42 


1 .470 


1.370 


- C . 1 0 C 


6.44 


1.030 


1.180 


0.150 



1.5 Amp Discharge 



Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 




0.00 


2.052 


2.064 


0.312 


0.7 5 


2 . 34 0* 


2 . 0 58 


0.31 8 


CD 

LO 

• 


2.020 


2 • C 4 9 


C . 3 29 


2.25 


2.000 


P .^36 


Q • C 7 6 


2.63 


1 . 99 G 


2.027 


C • n 37 


3.00 


1.980 


2.016 


C . 0 7 6 


3.38 


1.965 


2 . 0 0 G 


2 . 0 7 7 


4.13 


1 . Q 2 0 


1 .9 5? 


0.91 ° 


4.50 


1 . 8 8 0 


1.865 


_ p n •*, c. 

• t Aw 


4.65 


1.850 


1.813 


-0.037 


4.80 


1 . 8 1 0 


1 . 729 


-*'•081 


4.88 


1.770 


1.664 


- C . 1 ? 6 


4.95 


1 .720 


1.573 


-0.147 


5.0 3 


1.600 


1 . ^32 


_ r i 'p 


5.10 


1 . 1 0 c 


1.18 6 


0 . C ^ 6 


5.13 


0.850 


1 .058 


0 . 7 0 p 
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Numeric data corresponding to Figures k-25 and U-2 6 . 



COEFFICIENTS FOR THE MODEL 
SEE TABLE 4-3. 



3.6 Amp Discharge 




5.4 Amp Discharge 




Ah 


Volt 


Volt 


Diff 


Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 






(Exp) 


(Eqn) 




0 .0 G 


2 . C 17 


2.054 


0 . 2 37 


0.00 


2.014 


2.065 


0 . 0 5 1 


0.3 0 


2.010 


2.033 


0.023 


0.45 


1 .990 


^ ,4 0 1 

CL • ' ! 1 


c . r T i 


0.60 


2 . 0 C 0 


2 . C 12 


0.0 12 


0.90 


1 .970 


1.577 


0.007 


0.90 


1.984 


1 .991 


0.0 07 


1.35 


1.950 


1.933 


- 0.017 


1.2C 


1.970 


1.970 


C . 0 0 0 


1.8 0 


1.920 


1 .88 8 


- C .032 


1.50 


1.960 


1.948 


- 0 . 01 ? 


2.25 


1 . 89 0 


1.942 


- 0.0 4 8 


1.80 


1.945 


1 .926 


- 0 . 01 ° 


2. 70 


1.850 


1.794 


- n • r; 5 6 


2.10 


1.930 


1 .904 


- 0.026 


3.15 


1.800 


1 . 7^1 


- 0 . 0-9 


2.40 


1.910 


1.881 


- 0 . 0 ? c 


3.4 2 


1.750 


1.701 


-0 . r 49 


2.70 


1.890 


1.857 


- 0.033 


3.6 0 


1.700 


1.663 


- 0 . 37 


3.00 


1.860 


1.832 


- 0.028 


3.6 9 


1.660 


1.635 


- 0 . n 2 5 


3.3 0 


1.830 


1 .804 


- 0.026 


3.73 


1.570 


1 • 5 a l 


0 ♦ G 6 ’ 1 


3.60 


1.790 


1.769 


- 0.021 


3.87 


1 . 23 C 


1.4 = 0 


0 . 2 6 C 


3.78 ‘ 


1 .760 


1.741 


- 0.0 19 


3.96 


0.800 


0 . 7 50 


- 0 . r 5 0 


3.96 


1.710 


1.700 


- 0.010 










4.08 


1.650 


1.652 


o • n 02 










4.14 


1.590 


1.611 


0.021 










4.2 0 


1.480 


1 .544 


0 . 06 * 










4.26 


1 .280 


1 .402 


0.122 










4 . 3 ? 


0 .940 


0.8 75 


- 0.085 
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Numeric data corresponding to Figures 4-28 and 4-29- 



COEFFICIENTS FOR THE MODEL 



E S 


K 


C 


n 


Corr 


Ha 


Rb 


2 .020 


0.00128 


5-803 


1. 22270 


3.05102 


0 . 02100 


0 . 00022 



0.6 Amp Discharge 



Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 




0*00 


2.110 


2.0 19 


— 0 • Q 0 1 


1.80 


2.060 


1 .996 


-0.064 


3.60 


2.010 


1.973 


-0.037 


4.20 


1.990 


, 1 .965 


-0.025 


4.80 


1.970 


1 .957 


-0.013 


5.40 


1.910 


1.948 


0.C38 


5.70 


1.08 0 


1 . Q 42 


0.06? 


6.0 0 


1.830 


1 .935 


0.105 


6.12 


1.810 


1 .930 


0.120 


6.24 


1.770 


1 .923 


0. 153 


6.30 


1.740 


1 . a 1 6 


0.176 


6.36 


1.69 0 


1 .905 


0.215 


6.39 


1.63 0 


1.895 


0 . ? 65 


6.42 


1.470 


1 .873 


0 . 4 Q8 


6.44 


1.030 


1.359 


0.8 2° 



1.5 Amp Discharge 



Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 




0 . 2 0 


2.052 


2.01S 


- ? , n 1 4 


0.75 


2 . 0 4 C 


1 . c 94 


-C . C.46 


1*50 


2.020 


1 . ° 70 


-0.0^0 


2.25 


2.00 0 


1 . 9 4 6 


- c . r 5 4 


2.63 


1 . 990 


1.934 


-0.0 66 


3. OD 


1 .980 


1 .921 


_ n ore 
. ' . 


3 ♦ 3 8 


V. 963 


i .90S 


-0 . 055 


4*13 


1.920 


1.881 


- 0 . n 7 9 


4*50 


1.880 


1 •866 


- C • r 1 4 


4.65 


1 .850 


1.859 


0.038 


4.80 


1.810 


1 .84« 


0.039 


4 • 6 8 


1 . 770 


1 . 8 * 3 


0 . "• 73 


4.95 


1.720 


1 .835 


C.115 


5.03 


1.600 


1.825 


0 . 7 25 


5.10 


1.100 


1.839 


0.739 


5.13 


C .850 


1 . 8 : 1 


C . 9 5 1 
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Numeric Data corresponding to Figures — 3 0 and 1*-31. 
COEFFICIENTS FOR THE MODEL 

E s K C n Corr R a R b 

2.020 0.00128 5-803 1.222J0 3.05102 0.02100 0.00022 



3.6 Amp Discharge 



Ah 


Volt 


Volt 




(Exp) 


(Eqn) 


0.00 


2.017 


2.015 


0.3 0 


2.010 


1 .992 


0.60 


2.000 


... 1 .969 


0.90 


1 .984 


1.946 


1.20 


1 .970 


1.923 


1.50 


1.960 


1 *899 


1.R0 


1.945 


1 .876 


2.10 


1.930 


1 .852 


2.40 


1 .910 


1 .823 


2.70 


1.890 


1 .803 


3. GO 


1.860 


1 .778 


3.3 0 


1.830 


1 .751 


3.60 


1 .790 


1.721 


3.78 


1.760 


1 .699 


3.96 


1.710 


1 .670 


4.08 


1 • 65 0 


1 .640 


4.14 


1.590 


1.616 


4.20 


1 .480 


1.578 


4.26 


1 .280 


1.502 


4.3 2 


0 .94 Q 


1 .222 





5.4 


Diff 


Ah 



— 


— — 


- 0 . 0 02 


0 . 0 0 


- 0 . n 16 


a. 45 


- 0.031 


C .90 


- 0 . C 33 


1 .35 


1 

0 

. 

0 

& 


1.80 


- 0.061 


2.25 


- 0. 069 


2.7 0 


- 0.0 73 


3.15 


- 0 • 082 


3.4 2 


- 0 . C 87 


5.60 


- 0.062 


3.6 9 


- 0. 0 79 


3 • 7 8 


- C . 0 6 ^ 


3.8 7 


- 0.061 


3.96 


- 0.040 




- 0.010 




0.026 




0.098 




0.222 




0 .2 62 





Discharge 



Volt 


Volt 


Diff 


(Exp) 


(Eqn) 






2.014 


2 . C 1 2 


-0 . 


c 0 ? 


1.990 


1.963 


-0 . 


P 


1 .970 


1 .908 


-0 . 


0 62 


1.950 


1.856 


-0 . 


2 ° 4 


1.920 


1 . 3 0 2 


-0. 


! 1 0 


1.890 


1 .748 


_ 0 

• 


1 h 2 


1.850 


1 *692 


-0. 


1 53 


1 . 8 0 C 


1 .429 


-0 . 


1 71 


1.750 


1.5 3 3 


p 

• 


1 ^7 


1.700 


1.540 


— 0 . 


1 6 0 


1.660 


1 . 50S 


— 0 * 


1 52 


1.570 


I .457 


-n # 


115 


1.230 


1 . v ‘43 


f . 


1 15 



0.80 0 0 . 510 - 0 • ? ^ C 



Cm V 



7 0 s 



Numeric data corresponding to figures 4-33 and 4-34 



COEFFICIENTS FOR THE MODEL 

E g K C n Corr R a 

2.023 0.00771 5-803 I.2227O 1.39146 0.01540 0.00361 

R a R b 

0.01540 0.00361 



0.6 Amp Discharge 



Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 




0 .CO 


2.11C 


2.019 


-0 . C°7 


1 . s 0 


2.050 


1 .094 


-0 . r 66 


3 .60 


2.010 


r.979 


-0.040 


4 .20 


1.490 


1.9 6 0 


•0 .0 30 


4.80 


1 .970 


1.947 


-0.023 


5.40 


1.910 


-1 .9 2 C 


0 .015 


5 .70 


1 .8*0 


1.906 


0 .026 


6 .00 


1.830 


1 . ? 6 6 


0.036 


6.12 


1.910 


1.333 


0 .023 


6 .24 


1.770 


1.772 


0 .002 


6 .30 


1.740 


1.714 


-0.026 


6.36 


1 .690 


1.509 


-0 .CPI 


6.3° 


1.63C 


1.514 


-0.116 


6 .42 


1 .470 


1.350 


-0.120 


6.44 


1 .030 


1.15 3 


0.123 



1.5 Amp Discharge 



Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 




0.00 


2.052 


1 <“■ 
1 r 
1 • 
1 O 


-0 . 342 


0 .75 


2.04 0 


1.791 


-0 . C49 


1 .50 


• 'J £ 


1 ,°7 n 


-0 . 


2.25 


/ m J v- 


1.752 


— C . C - 9 


2 .*3 


1 .99 ' 


1.74? 


-0 .0 


3.00 


1.9? : 


1.731 


-0 . 147 


3.33 


1.963 


1.= 1° 


-0.7-15 


4.13 


1.920 


1 . °9< 5 


-0.132 


4.50 


1 . 8 f 0 


1.353 


-0.017 


4.63 


1.850 


1.347 


-0 .003 


4.80 


1.81' 


1 . 8 2 c 


0.015 


4.8°. 


1 .7 7 0 


1 ,°0 9 


0 .0 39 


4.9 5 


1.720 


1.7?7 


0.0*7 


5 .03 


1 . 600 


1 . ^54 


0 . 154 


5.10 


1 . 1 c 0 


I.O?" 1 


0 . *07 


5.13 


0.85 C- 


1 ,56« 


0.81° 
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Numeric Data corresponding to Figures 4-35 and 4-36 



Coefficients for the Model 



n 



Corr 



R> 



2.023 0.00771 5-803 1.22270 1.39146 0.01540 0.00361 



3.6 Amp Discharge 5.4 Amp Discharge 



Ah 


Volt 


Volt 


Diff 


Ah 


Volt 


Volt 


Diff 




(Exp) 


(Eqn) 






(Exp) 


(Eqn) 




0.00 


2.017 


7.902 


-0 .015 


9 .CC 


2.014 


1.906 


- 0 . r 15 


0 .30 


2.010 


1.385 


-0 .C25 


0.45 


1 .“9 c n 


^ .OS' 7 


-O .373 


0 .60 


O.OQC 


1.96 = 


-0.032 


0.50 


1.970 


1.9 19 


-0 . 0 c 1 


0.90 


1.984 


1 .950 


-9.034 


1 .35 


1.950 


1 .330 


-9.070 


1 .20 


1.970 


1.933 


-0.037 


1.80 


1 .9 29 


1 .049 


-0 .0°C 


1 .50 


1 .960 


1.915 


-0.045 


2.25 


1.899 


1.797 


-9.091 


1 .BO 


1.945 


1.897 


— 0 .043 


2 .79 


1.859 


1 .759 


- 0 .995 


2.10 


1 .930 


1.379 


-0 .0 51 


3.15 


1 . 8 c r, 


1.70 6 


-0.075 


2.90 


1.910 


1. 5 60 


-0.959 


3.42 


1 . 7 c 3 


1.565 


-0 . 755 


2.70 


1.890 


1.84 9 


-0.0B9 


3 .6 n 


1 .70 0 


1.625 


-0 .075 


3.00 


1.860 


1.819 


-0.941 


3.69 


1 .66 3 


1.593 


-0 . 9 67 


3.30 


1.830 


1 . 79 F 


-0.035 


= .■*£, 


1.57.9 


1.549 


-0 .0 3 9 


3.60 


1.790 


1.766 


-0.024 


3 .?•» 


1.23" 


1.417 


0.1 = 7 


3 .73 


1 .760 


1.743 


-0.C17 ’ 


3 .96 


9.39" 


'.492 


-O. 30 ? 


3 .9 6 


1.710 


1.703 


-o.ro 3 










4 .03 


1.650 


1.565 


0.015 










4.14 


1.590 


1.529 


0.039 










4 .20 


1.469 


1.579 


0.099 










4 .26 


1.2 90 


1 .44* 


0. 165 










4.32 


9.940 


0.982 


0 .042 
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